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We revisit the problem of electrons on a square lattice below half filling close to an Ising-nematic 
quantum critical point. For Fermi surfaces with sufficiently strong antinodal nesting, the static 
nematic susceptibility is maximal at the antinodal nesting wave vector within a simple RPA calcu¬ 
lation. We present a detailed analysis of the nematic susceptibility within Eliashberg theory and 
show that the strong interaction between Fermions in the antinodal regions shifts the maximum of 
the nematic susceptibility to slightly larger wave vectors. The corresponding order is akin to the 
incommensurate charge density wave with d-wave form factor found recently in some underdoped 
cuprate materials. At sufficiently high temperatures around T/t ~ 0.1 nematic fluctuations are 
strongest at zero wave vector. 


I. INTRODUCTION 

Incommensurate charge density wave (CDW) order has 
been observed recently in several cuprate high-T c super¬ 
conductors at small hole doping. 1- ' The measured or¬ 
dering wave vectors Q = (±Q o ,0) and (0, ±Q 0 ) are ori¬ 
ented along the principal axes of the Cu02 planes and 
take values around Q o ~ 0.25... 0.32 in reciprocal lat¬ 
tice units, corresponding to a modulation of the CDW 
order parameter with a period between three and four 
lattice spacings. Three recent experiments by Comin et 
al. [8], Fujita et al. [9] and Achkar et al. [10] showed 
that this CDW order has a predominant d-wave form 
factor in several cuprate families such as YBa 2 Cu306+ y , 
Bi 2 Sr2CaCu208+x and Ca 2 - x Na a ;Cu02Cl2. 

This observation of CDW order has prompted theo¬ 
retical interest in the origin of the axial ordering wave 
vectors Q = (±<2o,0) and (0, ±<9o)- Interestingly, many 
theoretical approaches predict a dominant instability 
towards a CDW with d-wave form factor at diagonal 
wave vectors (±QoiiQo)) whereas the experimentally 
observed axial wave vectors appear only as a subleading 
instability. 11-18 Other theoretical works argue that the 
dominant instability is indeed at the axial ordering wave 
vector, either because of the interplay between CDW 
and superconductivity, 19,20 the effect of a sufficiently 
large nearest-neighbor Coulomb repulsion V in a t-J-V 
model, 21 or because of other microscopic details. 22-25 

The observed charge-density wave with a d-wave form 
factor is akin to an incommensurate Ising-nematic or¬ 
der. For this reason we revisit the problem of two- 
dimensional metals close to an Ising-nematic quantum 
critical point (QCP), which has been studied thoroughly 
in the past, 26-3 ' with the aim to determine the wave 
vector at which nematic fluctuations are strongest. In 
fact, a similar analysis has been performed earlier within 
a random phase approximation (RPA) by Holder and 
Metzner, 12 where it was shown that the static nematic 
susceptibility is largest along 2 k F lines, i.e. at mo¬ 
menta connecting points on the Fermi surface with par¬ 
allel tangents. More precisely, they found that the ne¬ 


matic susceptibility for electron densities larger than 
van-Hove filling 38 is maximal at diagonal wave vectors 
Q = (±Q 0 > ±Qo) where two 2k F lines cross. Again, these 
results suggest that the ordering along axial wave vectors 
is a subleading instability. 

Here we demonstrate that the static nematic suscep¬ 
tibility can also be maximal at axial wave vectors Q = 
(±Qo, 0) an d (0, ±Qo); depending on the nesting proper¬ 
ties (i.e. the curvature) of the Fermi surface in the antin¬ 
odal region. Indeed, for sufficiently strong nesting of the 
Fermi surface in the antinodal region, the maximum of 
the nematic susceptibility is precisely at the antinodal 
nesting wave vector within a simple RPA calculation. 

Even though the antinodal nesting wave vector is closer 
to the experimentally observed ordering wave vector than 
the diagonal one, it misses one essential feature. Experi¬ 
ments suggest that the CDW ordering wave vector con¬ 
nects the tips of the Fermi arcs in the pseudo-gap regime, 
rather than the putative Fermi surface in the antinodes. 39 
This observation might be suggestive of the fact that the 
CDW order is a secondary instability of a genuine pseudo¬ 
gap phase which already features a reconstructed Fermi 
surface. 40,41 

Here we point out one possible resolution how the ob¬ 
served ordering wave vector could be reconciled with the 
picture of a CDW instability of an ordinary Fermi liquid 
with a large Fermi surface. In the nematic QCP scenario 
electrons close to the antinodes interact strongly and are 
no well-defined quasiparticles. Nematic fluctuations are 
strongest at 2 k F lines where quasiparticle coherence is 
important, however. It is thus conceivable that the max¬ 
imum of the nematic susceptibility shifts to larger wave 
vectors which connect points on the Fermi surface away 
from the antinodes where the quasiparticle coherence is 
larger, even though these points have no perfectly par¬ 
allel tangents. In this paper we present results from an 
Eliashberg-type calculation of the nematic susceptibility 
which suggest that this is indeed the case, even though 
this effect is rather small. The wave vector is only slightly 
shifted away from the antinodes, and we conclude that 
the nematic QCP scenario alone can hardly account for 


in Ref. [48], which is defined by the euclidean action 



FIG. 1: (Color online) Ordering wave vectors and Fermi sur¬ 
face. The black arrows indicate the antinodal nesting wave 
vectors, where the static nematic susceptibility is maximal 
within RPA for Fermi surfaces with sufficiently strong antin¬ 
odal nesting. In Eliashberg theory the maximum of the sus¬ 
ceptibility is shifted to slightly larger wave vectors away from 
the antinodes, indicated by the red dashed arrows. 


the experimentally observed CDW ordering wave vector. 
Superconducting fluctuations, which are clearly impor¬ 
tant in the context of cuprates, likely increase this shift 
of the CDW ordering wave vector, but a detailed study 
of this effect is beyond the scope of this paper. 

Another interesting result from the Eliashberg com¬ 
putation is that the maximum of the nematic suscep¬ 
tibility shifts to zero wave vector at sufficiently high 
temperatures above T/t ~ 0.1. The incommensu¬ 
rate CDW fluctuations at low temperatures thus ap¬ 
pear as ordinary, intra-unit-cell nematic fluctuations at 
higher temperatures. Indeed, intra-unit-cell nematic- 
ity has been reported in cuprates below the pseudogap 
temperature. 42-46 However, due to the gradual onset of 
CDW order as a function of temperature it is not clear 
at present if these observations are related, or if a regime 
with intra-unit cell nematicity exists above the onset of 
CDW order. 47 

The rest of the paper is outlined as follows: in Sec. II 
we introduce the model of electrons on a square lattice 
coupled to Ising-nematic order parameter fluctuations. 
Sec. Ill discusses RPA results for the nematic suscepti¬ 
bility and in Sec. IV we present results from Eliashberg 
theory. 


II. MODEL 


We start from the model of electrons on a square lattice 
coupled to an Ising-nematic order parameter introduced 
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Here the real field <j>k denotes the Ising-nematic order pa¬ 
rameter and s defines the distance to the critical point. 
Electrons with Matsubara frequency ui n , momentum q 
and spin a are denoted by the fermionic fields c 9tT , 
where we use the shorthand notation q = (w n ,q) and 
= /3 -1 Ibz ^ 2 ?/(2 7r ) 2 an d A is the interaction 
strength. In the following we use a standard tight bind¬ 
ing dispersion for the electrons with up to third-nearest 
neighbor hopping on a square lattice, which takes the 
form 


£ q = — 2f(cos q x + cos q y ) — 4 t' cos q x cos q y 

—2t"(cos2q x + cos2q y ) — fj, , (2.2) 

where /x is the chemical potential (the lattice constant is 
set to unity throughout this paper and in the following 
we measure energy in units of t = 1). 


On a square lattice the Ising nematic order has d x 2_ y 2 
symmetry and correspondingly the form factor Vk.q hr 
the interaction term of Eq. (2.1) is typically chosen as 
hk.q ~ G? q = cos q x — cos q y . In this case the nematic or¬ 
der parameter at finite k is identical to a bond-density- 
wave order with d-wave form factor which has been dis¬ 
cussed in the context of CDW ordering in underdoped 
cuprates. 9,13 By contrast, in the model introduced in 
Ref. [48] the Ising field fa lives on the sites i rather than 
the bonds of the square lattice. In the context of hole- 
doped cuprates this site-centered Ising field can be viewed 
as a deformation of a Zhang-Rice singlet which breaks 90° 
rotation invariance. 49 The form factor Vk.q in Eq. (2.1) 
corresponding to such a site-centered field fa takes the 
slightly different form 


I'k.q — dq^k /2 + dq—k /2 
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and depends explicitly on the fluctuation wave vector 
k, but obviously reduces to the standard definition at 
k = 0. In the following we will always use Eq. (2.3) in 
our computations. Note that this form factor respects 
all symmetries of an Ising-nematic order parameter with 
d x 2_ y 2 symmetry by construction, i.e. it is odd under 
7r/2-rotations and reflections about the nodal directions, 
and it is even under time-reversal as well as reflections 
about the x- and y-axis. Together with the Ising sym¬ 
metry (f>k —> —<j>k the action (2.1) is invariant under the 
square lattice point group. At the level of the bare ac¬ 
tion the Ising symmetry is broken for s < 0, in which 
case the point group symmetry is reduced from tetrag- 
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onal to orthorhombic. At a fixed ordering wave vector 
k the form factor (2.3) can be viewed as a d-wave order 
parameter admixed with an extended s-wave component 
s' = cos q x + cos q y . 

We emphasize that at the level of the bare action or¬ 
der parameter fluctuations are not momentum- and fre¬ 
quency dependent. Momentum- and frequency depen¬ 
dent terms will be generated by quantum fluctuations, 
however, which in turn determine the ordering wave vec¬ 
tor in the symmetry broken phase. 


III. NEMATIC SUSCEPTIBILITY IN RPA 

Within RPA the propagator Do(k, fi„) of the 0-field, 
which equals the nematic susceptibility, takes the form 

x (k, Q n ) = s- 2A 2 n 0 (k,n„) , (3.1) 

where the factor 2 arises from the electron spin. The wave 
vector at which the static susceptibility Ho(k, —>• 0 + ) 
is maximal determines the modulation wave vector of the 
nematic order parameter in the ordered phase. IIo(k, Q n ) 
is the bare d-wave polarization function defined by 

n 0 (k, f l n ) = — Go(k + q, Q n + w n )Go(q, co n ) V^ q+ k 

(3.2) 

where G'o(k, ui n ) = (iui n — £ k) 1 denotes the bare electron 
Green’s function. 

Holder and Metzner [12] computed a similar static d- 
wave polarization function for electrons on a square lat¬ 
tice using Vk.q ~ dq rather than Eq. (2.3) and showed 
that the nematic susceptibility is enhanced for momenta 
k satisfying the condition £(k+G)/2 = Oj with G a recipro¬ 
cal lattice vector. These special momenta connect points 
on the Fermi surface with parallel tangents and are simple 
lattice generalizations of 2/c^-momenta in systems with a 
circular Fermi surface. 50 At electron densities higher than 
van-Hove filling the maximum of the nematic susceptibil¬ 
ity is at the diagonal wave vectors Q = (±Qq F , ±Qq F ), 
where the two 2 kp lines with G = (2w,0) and (0, 27 t) 
intersect. A similar situation is shown in the right panel 
of Fig. 2, where we plot the static d-wave polarization 
function as defined above using Eq. (2.3) as a function of 
momenta at a temperature T/t = 0.01. Interestingly, 
these diagonal wave vectors connect antiferromagnetic 
hotspots on the Fermi surface. Indeed, after shifting mo¬ 
menta by (0, 27r) the condition of the two intersecting 2 kp 
lines £q /2 = £q/ 2 +tt = 0 is equivalent to the condition 
which determines the position of the antiferromagnetic 
hotspots £ qAF = ^qAp+Tr = 0, i.e. Q = 2(\af connects an¬ 
tiferromagnetic hotspots on opposite sides of the Fermi 
surface with parallel tangents. 

Here we want to point out that the momentum along 
the 2 kp lines where the nematic susceptibility is maxi¬ 
mal also depends on the curvature of the Fermi surface 
in the antinodal region around k = (0, n) and symmetry 
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FIG. 2: (Color online) Static d-wave polarization function 
IIo(k, 0) for different Fermi surface geometries. Left column: 
t!/t = —0.32, t"/t = 0.128, n/t = —0.85 (8% hole doping 
relative to half filling), right column: t'/t = —0.2, t"/t = 
0.05, n = —0.75 (12% hole doping). The upper graphs show 
density plots of no(k, 0) in one quadrant of the Brillouin zone 
and the insets indicate the shape of the Fermi surface. The 
lower graphs show cuts of no (k, 0) along the x-axis (blue thick 
line) and along the Brillouin zone diagonal (thin red line). 
For sufficiently strong antinodal nesting the maxima of the 
polarization function are at the axial wave vectors of the form 
k = (±Qo,0) and (0,±Qo), otherwise the maxima are at 
diagonal wave vectors k = (Qo,Qo) and symmetry related 
points. The grey dashed and dotted lines in the lower graphs 
show the d-wave polarization function as computed by Holder 
and Metzner [12] for comparison, i.e. using 14,q = 2d q instead 
of Eq. (2.3). 


related points in the model with a site-centered Ising- 
nematic field considered here. For the case of sufficiently 
strong antinodal nesting, the maximum is at a wave vec¬ 
tor Qq = (Qo.O) and symmetry related points, which 
connects opposite antinodal points of the Fermi surface, 
as indicated in Fig. 1. We show an example of this situ¬ 
ation in the left panel of Fig. 2, where we plot the static 
d-wave polarization function at ui = 0 and T/t = 0.01 
for a different Fermi surface geometry with less curva¬ 
ture in the antinodal region. The cuts along the prin¬ 
cipal axis and the diagonal of the Brillouin zone show 
clearly that the maximum is at Qo = (Q 0 , 0) rather than 
Q = (<2o F ,<2o F )- We emphasize that the form factor 
(2.3) which depends explicitly on the fluctuation wave 
vector k is crucial in order to observe this dependence on 
the Fermi surface geometry. For 14 q ~ e? q the maximum 
remains at the diagonal wave vectors in RPA, indepen- 
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FIG. 3: (Color online) Density plots of the static d-wave po¬ 
larization function IIo(k, 0) at temperature T/t = 0.01 (upper 
left) and T/t = 0.1 (upper right). The lower graph shows cuts 
along the x-axis and the Brillouin zone diagonal. Parameters: 
t'/t = —0.3, t"/t = 0.1, n/t = —0.8 corresponding to 6.6% 
hole doping. 


dent of geometric details of the Fermi surface. 

At finite temperature the 2 Icf singularities in the ne¬ 
matic susceptibility are broadened and the maximum is 
at a slightly smaller wave vector, as shown in Fig. 3, 
where we compare the static d-wave polarization func¬ 
tion at T/t = 0.01 and T/t = 0.1. We note that this is 
no longer the case if effects beyond RPA are taken into 
account. Within Eliashberg theory the maximum shifts 
to Q = 0 at T/t = 0.1, as will be shown in Sec. IV. 

As mentioned in the introduction, the experimentally 
observed CDW ordering wave vectors do not seem to con¬ 
nect the Fermi surface at the antinodal points. Rather, 
it has been suggested that the ordering wave vectors con¬ 
nect the tip of the Fermi arcs in the pseduogap phase. 39 
Even though the RPA results above give the correct qual¬ 
itative form of the ordering wave vectors if the Fermi sur¬ 
face is sufficiently nested in the antinodes, there is clearly 
something missing. Indeed, nematic fluctuations give rise 
to a strong interaction between electrons in the antin¬ 
odes, which in turn leads to a suppressed quasiparticle 
coherence for electrons with momenta in the antinodal 
region. The nematic susceptibility is very sensitive to 
quasiparticle coherence, however, and thus it is reason¬ 
able to expect that the maximum of the susceptibility 



FIG. 4: (Color online) Static d-wave polarization function 
IIo(k, 0) at temperature T/t = 0.01 along the k x axis ( k y = 0). 
The thick solid line shows the result using the renormal¬ 
ized electron propagators (3.3) with the quasiparticle residue 
Zk from Eq. (3.4) and has been multiplied by a factor of 
5 for better visibility. The thin dashed line shows the re¬ 
sult with Zk — 1 for comparison. The reduced quasiparticle 
residue in the antinodal region shifts the local maximum of 
the susceptibility from Qo = (1.31,0) to a larger wave-vector 
Q = (1.40,0), indicated by the red arrow in Fig. 1. Parame¬ 
ters: t'/t = —0.3, t"/t — 0.1, ^ = —0.8. 


shifts to larger wave vectors connecting points on the 
Fermi surface that are no longer perfectly nested, but 
have a larger quasiparticle residue. For this reason it is 
important to study the susceptibility beyond RPA, where 
such effects are not taken into account. The next section 
contains detailed results for the nematic susceptibility 
within Eliashberg theory, but for now we keep our discus¬ 
sion at a phenomenological level. The simplest approach 
to include electronic quasiparticle coherence in RPA is 
to compute the d-wave polarization function in Eq. (3.2) 
with renormalized electron propagators 

G(k,cj n )= . Zk c . (3.3) 

n sk 

The quasiparticle residue Z k is unity along the nodal 
lines k x = ±fc y , where the electrons do not interact with 
nematic fluctuations, and falls off towards the antinodes. 
A very simple phenomenological parametrization of the 
quasiparticle residue consistent with the square lattice 
symmetry takes the form 

Zk = 1 — (cos k x — cosky ) 2 /4 . (3.4) 

In Fig. 4 we show the static d-wave polarization function 
computed with a quasiparticle residue of this form. As 
anticipated, the local maximum of the nematic suscepti¬ 
bility indeed shifts to larger wave vectors, indicated by 
the red dashed arrows in Fig. 1. Including Z k the wave 
vector of the maximum changes from Q = (1.31,0) to 
Q = (1.40,0), which connects points on the Fermi sur¬ 
face away from the antinode with tangents that are no 
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longer perfectly parallel. 

It is important to note that the global maximum of 
the susceptibility changes from axial to diagonal wave 
vectors Q = (Qq F ,Qq F ) in this simple scenario consid¬ 
ered here. This is not surprising, because the quasipar¬ 
ticle residue at the antiferromagnetic hot-spots is larger 
in the ansatz (3.4) and the diagonal wave vector still 
connects points on the Fermi surface with parallel tan¬ 
gents. Within Eliashberg theory, where the damping of 
electronic quasiparticles along the Fermi surface is deter¬ 
mined selfconsistently, this is not necessarily true, as will 
be shown in the next section. 


IV. ELIASHBERG THEORY 

We now turn to a self-consistent one loop calculation of 
the d-wave polarization function II(k, fl n ) and the elec¬ 
tronic self-energy £(k,w n ). In the spirit of Eliashberg 
theory vertex corrections are neglected and we comment 
on the validity of this approximation later. Writing the 
full electron and nematic fluctuation propagators as 

G-\k) = iw„-&-E(fc) (4.1) 

D-^k) = s - 2A 2 II(fc) , (4.2) 

the coupled Eliashberg equations for the electronic and 
bosonic self-energies take the form 

£(fc) = A 2 ^G(fc-g)£( 9 )V q 2 k _ q /2 (4.3) 

Q 

H(fc) = -Y,G(k + q)G(q)Vl q+k/2 ! (4-4) 

where we used the abbreviation k = (w n ,k) again. We 
solve these coupled equations numerically after analytic 
continuation to real frequencies by a fixed point itera¬ 
tion on a grid of up to 210 x 210 points in the first 
Brillouin zone and 201 points along the frequency axis. 
For all results shown in this section the coupling con¬ 
stant is chosen as X/t = -/(hi ~ 0.316 and the gap 
s = s — 2A 2 max k ReII(k, 0) of the nematic fluctuation 
propagator D(k, w) is fixed at the value s/t = 0.01 at a 
temperature T/t = 0.01 (we do not fix the gap at s/t = 0 
at T = 0 for numerical reasons). At higher tempera¬ 
tures the boson gap is determined self-consistently and at 
T/t = 0.1 we obtain s/t = 0.022. Furthermore, we fixed 
the renormalized chemical potential jl = /r — Re T*(kp, 0) 
at the nodal point to avoid large shifts of the Fermi sur¬ 
face. 


A. Nematic susceptibility 

In Fig. 5 we plot the static d-wave polarization func¬ 
tion II(k,u; ra —> 0+) for two different temperatures as a 
function of momenta. In contrast to the RPA results 
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FIG. 5: (Color online) Density plots of the static d-wave po¬ 
larization function Ft(k, 0) computed within Eliashberg the¬ 
ory at a temperature T/t = 0.01 (upper left) and T/t = 0.1 
(upper right). The lower panel shows cuts through II(k, 0) 
along the x-axis and the Brillouin zone diagonal. Parameters: 
t'/t = -0.3, t"/t = 0.1, ji/t = -0.85. 


a substantial broadening of the 2k f singularities is al¬ 
ready visible at low temperatures T/t = 0.01, which we 
attribute to the the reduced quasiparticle coherence of 
electrons close to the antinodes, as discussed earlier. The 
maximum of the static nematic susceptibility remains at 
the axial wave vectors Q = (1.20, 0) and symmetry re¬ 
lated points. We note, however, that it is almost degener¬ 
ate with the local maxima at the diagonal wave vectors, 
where the susceptibility is a quarter percent smaller. 

The wave vector Q = (1.20,0) is slightly larger than 
the antinodal nesting wave vector and connects points 
on the Fermi surface away from the antinodes. In¬ 
deed, we determined the antinodal nesting wave vec¬ 
tor from the maximum of the electron spectral function 
A(k,w) = — 4lmG(k, ico n — > uj + *0+) at T/t = 0.01 and 
zero frequency to be Q 0 = (1.18, 0) and the wave vector 
Q thus connects hotspots at k = (±0.60, 2.79) on the 
Fermi surface. 

Interestingly, the 2 kp singularities are no longer vis¬ 
ible at higher temperatures above T/t ~ 0.1 and the 
nematic susceptibility is maximal at k = 0 instead. The 
appearance of such a preemptive k = 0 nematic order 
has been discussed earlier in the context cuprates and 
iron-pnict ides. 1 9 ’ 1 
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FIG. 6: (Color online) Electron spectral function A(k, to) at a 
temperature T/t = 0.1, plotted as a function of frequency for 
three different momenta along the Fermi surface. The inset 
indicates the position of the three momenta on the Fermi 
surface. Parameters are the same as in Fig. 5. 


B. Electron spectral function 

Finally we present results for the electron spectral 
function A(k, w) within Eliashberg theory. In Fig. 6 we 
plot A(k, w) at a temperature T/t = 0.1 as function of 
frequency for three different momenta along the Fermi 
surface. The spectral function exhibits a peak at the 
Fermi energy and its width increases towards the antin¬ 
odes. There is clearly no pseudogap, i.e. no suppression 
of spectral weight at the Fermi energy close to the antin¬ 
odes. Similar results were obtained already in Ref. [36]. 
It has been argued, however, that the interplay between 
CDW and d-wave superconducting fluctuations leads to 
the formation of a low-temperature pseudogap. 52 


C. Validity of Eliashberg theory 

The validity of the Eliashberg approximation for two- 
dimensional metals close to a k = 0 nematic critical point 
has been analysed in detail in Ref. [33]. One of the re¬ 
quirements there is that the effective four-fermion inter¬ 
action A 2 I1(Q, 0) is small compared to the bandwidth, 
i.e. A 2 /s <C t in order for the leading order vertex cor¬ 
rection shown in Fig. 7 to be small. Here we want to 
emphasize one important difference for the case of ne¬ 
matic fluctuations at the axial wave-vector Q = (Qo?0)- 
In fact, the scattering geometry in this case is such that 
one of the internal electron propagators is always off- 
shell and the first order vertex correction is small even 
in a regime where A 2 /s > t as considered here. Indeed, 
setting p = k = Q in the diagram in Fig. 7 one can 
easily see that the electron with momentum q + 3Q/2 
is far away from the Fermi surface if the incoming and 
outgoing electrons are at the hotspots connected by Q. 



FIG. 7: Lowest order vertex correction which is neglected in 
Eliashberg theory. Solid lines denote the electron propagator 
G(k), wiggly lines represent the nematic fluctuation propaga¬ 
tor D(k). 


To be more explicit, the vertex correction shown in 
Fig. 7 takes the form 

r K = X 3 ^2 G (p + q + k/2)G(p + q-k/2)D(p) 

V 

xVk.q+pFp q+ p=kF p q+E ±k . (4.5) 

Note that G(k) and D(k) denote the fully dressed elec¬ 
tron and fluctuation propagators obtained within Eliash¬ 
berg theory. In order to compute the renormalization of 
the coupling constant A we project the vertex correction 
-^kq zer0 external frequencies onto the d-wave form 
factor 14.q from Eq. (2.3) and obtain 

“-sEwi". <«> 

k.q 

where the factor 1/2 is for normalization. Computing 
the integral we obtain S A/A = —0.255 at T/t = 0.01 and 
SX/X = —0.318 at T/t = 0.1. Note that these vertex cor¬ 
rections are relatively small even though the four-fermion 
interaction A 2 /s is larger than the bandwidth t by a fac¬ 
tor of 10 for the parameters used here. 


V. CONCLUSIONS 

We presented numerical results for the nematic suscep¬ 
tibility of electrons on a square lattice below half filling. 
In contrast to models of electrons coupled to a bond- 
centered Ising-nematic order parameter field, the model 
used here exhibits a site-centered field, which can be in¬ 
terpreted as a deformation of a Zhang-Rice singlet cen¬ 
tered on the copper sites in the context of hole-doped 
cuprates. Within RPA the nematic susceptibility is max¬ 
imal either at diagonal-wave vectors connecting antiferro¬ 
magnetic hotspots on the Fermi surface, or at the antin- 
odal nesting wave vectors if the curvature of the Fermi 
surface at the antinodes is sufficiently small. In order 
to observe this dependence on the Fermi surface geom¬ 
etry it is crucial that the Ising-nematic order parameter 
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has a form factor with d x 2_ y 2 symmetry which depends 
explicitly on the fluctuation wave vector. We argued 
that the strong coupling between electrons and nematic 
fluctuations close to the antinodes reduces the electronic 
quasiparticle coherence and shifts the maximum of the 
nematic susceptibility to larger wave vectors at low tem¬ 
peratures. Results from an Eliashberg-type computation 
support this picture but indicate that this shift is rather 
small and can hardly account for the experimentally ob¬ 
served CDW ordering wave vectors in materials such as 
Bi2201, where the CDW ordering wave vector is roughly 
twice as large as the antinodal nesting wave vector. 8 In¬ 
cluding the effect of d-wave superconducting fluctuations 
potentially increases this shift of the ordering wave vec¬ 
tor, because such fluctuations reduce the electronic spec¬ 
tral weight at the Fermi level in the antinodal region. We 
leave this problem open for future study. 

Moreover, we showed that the nematic susceptibility is 
maximal at zero wave vector above temperatures on the 
order of T/t = 0.1. This seems to be in agreement with 
experiments on underdoped cuprates, where observations 
of intra-unit cell nematic order have been reported be¬ 
low the pseudogap temperature, which is parametrically 
larger than reported onset-temperatures of CDW order. 


It is important to emphasize, however, that the gradual 
onset of CDW order makes it difficult to determine such 
an onset-temperature and at present it is not clear if a 
regime with intra-unit cell nematic order exists above the 
onset of CDW order in the cuprates. 47 

Lastly, even though the wave vectors obtained in the 
Eliashberg calculation are in qualitative agreement with 
experiments on underdoped cuprates, the simple scenario 
of electrons close to a nematic quantum critical point 
clearly doesn’t describe the pseudogap behavior, as we 
do not observe a suppression of electronic spectral weight 
at the Fermi energy in the antinodal region. 
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